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Abstract. This work explores the possibility to transfer the parity law of the singlet-triplet gap established 
for square ladders (gapped for even number of legs, gapless for odd number of legs) to fused polyacenic 
1-D systems, i.e., graphite ribbons. Qualitative arguments are presented in favor of a gapped character 
when the number n u of benzene rings along the ribbon width is odd. A series of numerical calculations 
(quantitative mapping on spin 1/2 chains, renormalized excitonic treatments and Quantum Monte Carlo) 
confirm the parity law and the gapless character of the ribbon for even n u . 

PACS. 71.10.-W Theories and models of many-electron systems 71.15.Nc Total energy and cohesive energy 
calculations 75.10.Jm Quantized spin models 



1 Introduction 

In the recent past the properties of some quasi 1-D strongly 
correlated materials, namely cuprate ladders, have attracted 
much interest from solid state physicists |l|. As a major 
result they have established that the spin gap (i.e., low- 
est singlet to triplet excitation energy) of spin 1/2 ladders 
presents a parity law: the ladders are gapped (have a fi- 
nite excitation energy) when the number of legs is even, 
and gapless (degenerate singlet and triplet state) for odd 
number of legs 0. This result was not expected, since one 
may consider the ladders as intermediate between the sim- 
ple 1-D chain and the square 2-D lattice, which are both 
gapless. 

Finite ribbons of fused polyacenes (with CH bonds on 
the most external doubly-bonded carbons) can be seen as 
organic analogs of the cuprate ladders. They are quasi 1-D 
fragments of graphite, which is gapless. The non-dimerized 
linear polyene is also known to be gapless. In light of the 
recent investigations on both carbon and spin nanotubes, 
we can extend the discussion on the link between magnetic 
properties and geometry of such objects pj]. In particular, 
are the fused polybenzenoid ribbons always gapless ? 

Although the 7r-electrons of the conjugated hydrocar- 
bons are not strongly correlated, it has been shown twenty 
years ago that they can be treated accurately through 
S = 1/2 Heisenberg Hamiltonians 



H 



(i,3) 



(1) 



i.e., as spins interacting with their nearest neighbours 
through antiferromagnetic (AF) couplings J > 
This effective model is not based on the usual perturba- 
tive approach in the strong coupling limit; instead the 
amplitude of the exchange integral is given by the ex- 
act solution of the 2-site problem. More explicitely, for 
a Hubbard Hamiltonian, the AF coupling is no longer 
taken as its perturbative estimate J = 4t 2 /U, but rather 
as the exact energy differen ce between the triplet and 
singlet states J = -(U - VU 2 + 16i 2 )/2. This view of 
7r-electron systems offers simple rationalizations of many 
of their properties. In particular, a geometry-dependent 
Heisenberg Hamiltonian (with J(r) depending on the dis- 
tance r between sites) has been extracted from accurate 
calculations of the singlet and triplet states of ethylene 
and happens to be a quantitative tool for the ground and 
lowest excited states of conjugated hydrocarbons [HIB"] (a 
spin-independent V(r) potential is of course also needed 
in that case to take into account the r-dependences of the 
localized cr-bonds). This analysis was proved in various 
situations to be very predictive: for instance, the ground- 
state geometries of a large series of hydrocarbones are ac- 
curately reproduced by minimizing the lowest eigenvalue 
of the geometry-dependent Hamiltonian with respect to 
the bond distance |3j. It was also shown that the vertical 
excitation spectrum to the lowest triplet states is in good 
agreement with experiments [$] and that the excited state 
geometries and vertical emission energies agree with 
accurate ab initio calculations. This magnetic model has 
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Fig. 1. Two-leg ladder (A) and n u = 1 fused polyacene (B), 
and their mapping into dimerized spin 1 /2 chains. In this figure 
and the following ones, the symbols o and * correspond to 
up and down spins respectively; the resulting effective models 
appear just below each lattice. We denote with a single- and 
double- line the effective couplings when they are different. For 
the graphitic ribbons, dotted lines separate the different blocks 
used in the calculation (see text). 



been widely and successfully used by Robb et al. un- 
der the acronym MM-VB (Molecular Mechanics Valence 
Bond) for the study of the photochemistry of conjugated 
hydrocarbons. 

In fact, it is quite natural to expect, at least for bipar- 
tite lattices, the lowest excitation energies to be correctly 
described by a (properly normalized) Heisenberg Hamilto- 
nian as the lowest states are dominated by neutral VB con- 
figurations and as the antisymmetrization favours an anti- 
ferromagnetic order independently of the value of U • Of 
course, the charge excitations, essentially leading to ionic 
states, are not accessible with these approaches, but it is 
known for idealized ID chains [Hj and from experiments 
on conjugated molecules [Hj that the dipolarly-allowed VB 
ionic states are located at higher energies than the lowest 
triplet (and even than neutral singlet states). 

We now only consider the purely magnetic model of 
Eq. JQ. A few years ago density matrix renormalization 
group (DMRG) calculations have been reported, concern- 
ing the singlet-triplet gap of the simplest polyacenic chain, 
built of aligned fused benzene rings jT0|. The extrapolated 
calculated gap is finite and close to 0.1J. Actually the 
polyacene can be viewed as a two-leg ladder in which one 
rung over two has vanished (see Fig.^J. One might wonder 
whether there is a similarity between the three-leg ladder 
and a fused polyacenic infinite ribbon with two ranks of 
benzene rings (see Fig. El - From qualitative arguments one 
may conjecture that the singlet-triplet excitation in such 
polyacenic ribbons of graphite is finite when the number of 
superposed rings in the width of the ribbon, n w , is odd and 
vanishes when is even. Numerical calculations using the 
quantitative mapping on known 1-D chains, renormalized 
excitonic method (REM) pi] , and Quantum Monte Carlo 
calculations (QMC), show that the gap indeed vanishes 
for even n u and is non-zero for odd n u ribbons. 



2 Qualitative arguments 

Qualitative arguments can be used to rationalize the par- 
ity law of ladders, which can also be applied to 1-D fused 
polyacenes. They are based on the real-space renormaliza- 
tion group (RSRG), originally proposed by Wilson [12] . 
Concerning spin lattices one may consider that they are 
built from blocks rather than from sites. Of course these 
blocks interact. If the blocks do not have a singlet ground 
state, but a doublet, or a triplet, they can be seen as in- 
teracting effective spins ^3]. A quasi 1-D lattice can then 
be easily transformed Ql] into a simple 1-D spin chain, 
the properties of which are well known. Considering 2-leg 
ladders, one may define (2N a + 1) sites blocks, which have 
a ground state doublet (S z = ±1/2). But these blocks do 
not have equal interactions with their left and right near- 
est neighbors (cf. figure (A)). Consequently, the 2-leg 
ladder maps into a dimerized (i.e., bond-alternated) spin 
1/2 chain, which is known to be gapped. An alternative 
partition uses blocks with even number of sites (figure El 
(A)) which have a triplet ground-state. The resulting ef- 
fective S = 1 chain exhibits the famous Haldane gap ^S] 
and therefore both partitions predict a gap. 

If one applies the same arguments to the n u = 1 poly- 
acenic chain one obtains a similar mapping pictured in 
figureni(B). The partition into blocks of (4iVs + 1) sites 
and (47Vs + 3) sites produces an alternating dimerized 
spin chain, which is gapped. The partition into 4/Vs sites 
blocks with triplet ground states (figure leads to a Hal- 
dane gap Both partitions suggest a finite excitation 
energy. 

The ladder with odd number of legs can be partitioned 
into blocks with odd number of sites which have equal 
interactions with the left and right neighbors, and the 
ladders can be mapped into a non-dimerized gapless 1-D 
S = 1/2 chain, (cf. figure (A)). For the n u = 2 polyben- 
zenoid ribbon the simplest partition defines 9-site blocks 
presenting a S — 1/2 ground state and equal AF interac- 
tions with left and right nearest neighbors (cf. Fig. 0(B)). 
The resulting mapping to an effective AF uniform chain 
indicates that the = 2 polyacenic ribbon should not 
be gapped. This is actually one of the two possibilities 
allowed by the famous Lieb-Schultz-Mattis theorem jlfij 
which can be applied for even n w (where the unit cell 
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Fig. 2. Mapping of ladder and polyacene into a spin-1 chain. 
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Fig. 3. Three-leg ladder (A) and n u = 2 fused polyacene (B), 
and their mapping into spin 1/2 chains. 

contains an odd number of spin 1/2). If even n w lattices 
do have a gap, this theorem guarantees that the system 
breaks translation symmetry. 

The arguments can easily be generalized to any thick- 
ness of the ribbon according to figure which proposes 
a mapping into an uniform S — 1/2 spin chain for even 
values of n w and into an uniform integer-spin 5 chain for 
odd n^j. Of course the 5 — 1 chain may be transformed 
into a dimerized chain of 5 — 1/2 blocks by shifting one 
external carbon from one block to its right side neighbor 
block. 




(A) (B) 

Fig. 4. Generalization of the mappings on non-dimerized AF 
5=1/2 chains for even (A) and S = 1 chains for odd 
(B). Moving the (*) atom to the right-side block leads to an 
effective dimerized S = 1/2 AF chains. 

3 Numerical studies 

We now turn to numerical verifications of that conjecture 
for riu, — 1, 2 and 3 with several different techniques. Note 



that DMRG calculations have already been performed for 
the n u = 1 case dJJ. For larger n w , the number of states 
to be kept becomes too large and the computation cost 
is prohibitive. Therefore, after providing quantitative es- 
timates for the spin gap using approximate techniques, 
we have adressed this issues by using an efficient QMC 
method. 



3.1 Mapping into 1-D chains 

For riu) = 1 the chain may be considered as built from 
A and B blocks of 7 and 9, 9 and 11, 11 and 13, or 13 
and 15 sites, according to figuren](B). One then obtains a 
dimerized 5=1/2 spin chain with different values J\, J2 
between the A — B and B — A blocks. The effective inter- 
actions Ji and J2 are directly obtained from the energy 
difference between the lowest triplet and singlet energies 
of A - B and B - A blocks : 



Blocks 


7-9 


9-11 


11-13 


13-15 


Ji 


0.221 


0.196 


0.181 


0.170 


■h 


0.166 


0.123 


0.095 


0.075 



A previons study of the dimerized 5 = 1/2 chain has 
suggested that the gap follows the law : 

AE = (J 1 + J 2 )S - 71 with ,5=^— ^. (2) 

[J 1 + Ji) 

Applying Eq. Q one obtains 

AE = 0.097J for the (7-9) blocks, 
AE = 0.112J for the (9-11) blocks, 
AE = 0.121J for the (11-13) blocks, 
AE = 0.126J for the (13-15) blocks, 

hence a finite gap. 

For n u = 2 the simple mapping into a non-dimerized 
S = 1/2 AF spin chain of 9-site blocks (cf. figure (B)) 
pleads in favor of a gapless character, but this construction 
assumes only nearest-neighbor effective interaction Jnn, 
neglecting for instance next-nearest neighbors interactions 
Jnnn- In fact, an uniform AF chain becomes gapped 
when the ratio Jnnn / Jnn is larger than 0.241 jT7) . We 
have extracted the effective couplings between 9-site blocks 
(figure (B)) from the exact spectrum of the trimer of 
blocks and we found an AF coupling between NN blocks 
Jnn — 0.15896 (in good agreement with the value ex- 
tracted from the dimer Jnn — 0.16622) and a surprisingly 
large ferromagnetic J' NNN = —0.1705 coupling. Since an 
AF 5=1/2 chain with ferromagnetic coupling between 
NNN sites is not gapped, the n u = 2 ribbon should be 
gapless. 

For = 3 one may define a dimerized AF 5=1/2 
chain of 13 and 11 sites respectively as pictured in figure 
El One obtains two values of the inter block AF coupling 
(Ji = 0.115, J2 = 0.108). Using equation (J2J one obtains 
AE{n^ = 3) = 0.018J. 
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Fig. 5. Definition of the blocks for a mapping into a dimerized 
5=1/2 chain of n u = 3 ribbon. 



3.2 Renormalized excitonic calculations 

The recently proposed renormalized excitonic method 
is again based on a periodic partition into blocks, but now 
the blocks have an even number of sites, a singlet ground 
state and a triplet lowest excited state ij)* A . One defines 
a model space for the AB dimers made of two adjacent 
blocks, spanned by local singly excited states V'aV'b an d 
ipj^ip%- Knowing the spectrum of the AB dimer, it is pos- 
sible to define 

- the effective energy of V'aV'b an d iPa^Bi 

- the effective interaction between them. 

These informations allow to apply the excitonic method. 
For the infinite lattice, this method leads to the following 
expressions of the excitation energy 



AE°°(N S ) = 2AE(2N S ) - AE(N S ) 



(3) 



obtained from the gap values for a TVs-site block (AE(Ng)) 
and for a dimer (AE(2N S )). 

The method has been applied with success to the 2- 
leg ladders and can be applied to = 1, n w = 2 
and = 3 polybenzenoid ribbons using the design of the 
blocks pictured in figure El 

We find a finite gap for n w = 1 polyacene: 

AE°°(8) = 0.068 J for N s = 8 sites blocks, 
AE°°(12) = 0.094 J for N s = 12 sites blocks. 

Since the Ng 1 components of the excitation energy cancels 
in the expression of AE°°(Ns), an extrapolation in terms 
of A^ 2 leads to AE(n u = 1) = 0.103J. 

The result of the REM method for the n w = 2 lattice 
from the Ns — 12 sites blocks is one order of magnitude 
smaller, Z\_E°°(12) = 0.013 J. Extrapolation is not possible 
in this problem, but this result supports the conjecture 
that this fused polybenzenoid ribbon is not gapped. 

For — 3 unequal blocks of Ns = 14 and Ns = 10 
sites have been used, according to figureHU Since the blocks 
are different it is necessary to generalize the algebra of 
Ref. 6, as performed in the Appendix. The final result is 
AE°°(14, 10) = 0.0125 J. Even if the results for the gaps 
for = 2 and 3 are quite close, these calculations support 
the idea of a constrast between the odd and even ribbons. 




Fig. 6. Definition of blocks for the calculation of the energy 
gap through the REM. 



In order to confirm this conjecture, we now present accu- 
rate Quantum Monte Carlo results for the n w = 1,2 and 
3 ribbons. 



3.3 Quantum Monte Carlo 

As the spin lattices are not frustrated, efficient QMC al- 
gorithms are available which allow to simulate with high- 
precision large systems at finite, albeit extremely small, 
temperatures. Here we use a multi-cluster continuous time 
[Tfl] loop algorithm [2Qj which is free of any systematic er- 
rors. We simulate systems of sizes up to length L = 800 
(the total number of spins is N = 2(n w + 1)L) and at 
inverse temperature up to (3 J = J/T = 1000. We use pe- 
riodic boundary conditions along the legs. To determine 
the gapfull/gapless nature of the systems, we calculate the 
correlation length in imaginary time £ T with the help of a 
standard second moment estimator [21,22 



2tt 



X {UJ = 2TT/P) 



1/2 



(4) 



where x{ u ) = Jo dTe wr \(r) is the Fourier-transform of 
the imaginary time dynamical antiferromagnetic structure 
factor X (t) = ^ E^R^ jg dtS?(t)S/(t + r). It can 
be shown that for a gapped system, £ T converges to the 
inverse spin gap in the thermodynamic limit at zero tem- 
perature : linii ^oo £ T (L, /3) = (AE)^ 1 . On the other 
hand, when the system is gapless, S,^ 1 is an upper bound 
of the finite-size gap for any finite L and (3, In fig. we 
represent the inverse of the imaginary correlation length 
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Fig. 7. Quantum Monte Carlo results for the inverse of the 
imaginary time correlation length as a function of temperature 
for the n w = 1 (bottom panel), n„ = 2 (middle panel) and 
n-u, = 3 (top panel), for different system sizes L. Values of the 
extracted gaps AE(n w = 1) = 0.1112(7) J and AE(n u = 3) = 
0.0164(4) J are represented by dashed-lines for n w = 1 and 3. 



(^Ct)" 1 as a function of the temperature T/J in log-log 
scale, for the three types of ribbons n w = 1, 2 and 3. This 
representation is useful to see if the system is gapped as 
(J£t) saturates at low temperatures to the gap value 
AE/J. 

For flu = 1, (J^t)^ 1 clearly converges at low T to a 
minimum value identical for system sizes L — 100 and L = 
200, indicating that finite-size effects are absent. From the 
results for the largest L at the lowest T, we extract the 
value of the spin gap AEfo^ = 1) = 0.1112(7) J, in perfect 
agreement with DMRG calculations [10] . 

We find no saturation of (J£ r ) _1 at the lowest tem- 
perature for the largest system size L = 800 studied for 

= 2. Data for smaller systems (L = 100, 200) present 
signs of saturation at low enough T towards values which 
depend on the system size: this is naturally interpreted as 
the signature of finite-size gaps. Strictly speaking, the nu- 
merical data for the largest L at the lowest T — 0.001 can 
only put an upper bound AE(n u = 2) < 0.0030J on the 
value of the gap. However, the general form of the T de- 
pendence of (J£ T ) -1 and the fact that we still have finite- 
size effects for large systems at low T naturally indicate 
that the n w = 2 system is gapless, i.e. AEfa^ = 2) = 0. 

Finally for n w = 3, a convergence of (J£ r ) -1 is recov- 
ered at low enough T towards a size- independent constant. 
As for the n u = 1 case, this indicates that the system 
is gapped and we obtain the gap value AE(n u = 3) = 
0.0164(4) J. Please note that we had to resort to large sys- 
tems at very low temperatures (L = 800 and T = 0.001) 
to assert the complete convergence of our data: this is as- 
cribed to the small value of the gap (an order of magnitude 
lower than for n w = 1). 

In conclusion, the QMC simulations unambiguously 
proove the even/odd number n w effect for the gapless/gapfull 



nature of 1-D fused polyacenic ribbons. As for spin lad- 
ders, we find that the value of the gap decreases with the 
number n w of rings along the width for odd n w ribbons. 



4 Conclusion 

The conjecture of the existence of parity law concerning 
the spin gap in polybenzenoi'd ribbons (vanishing / finite 
gap for even / odd n^) was proposed on the ground of 
qualitative arguments. This conjecture is confirmed by a 
consistent set of numerical calculations summarized in ta- 
ble l.The REM fails to give a zero gap for n w = 2 but 
indicates a contrast between odd and even ribbons. The 
quantitative mapping and the accurate QMC calculations 
confirm the gapless character of the = 2 ribbons and 
agree on the amplitude of the gap for n u = 1 and n w = 3. 
For graphite ribbons [4,5J, for which J ~ 2.2 eV at the 
typical r cc distance (1.395A), the gaps should be close to 
0.23 eV for n u = 1 and 0.03 eV for n u — 3. Despite the 
semi-metallic nature of graphite, it seems unlikely that 
the odd-n w ribbons present a zero charge gap, below their 
finite spin-gap. We finally note, in analogy with the study 
of the influence of four-spin operators in spin ladders |23j, 
that it would be worth considering the possible impact 
of six-spin operators on such graphitic ribbons since their 
amplitude is not negligible in hexagons [24H25] . 





Mapping REM QMC 


n u = 1 


0.126 0.103 0.1112(7) 


n u = 2 


0.013 


n u = 3 


0.018 0.0125 0.0164(4) 



Table 1. Calculed gaps (in units of J) obtained from different 
methods. 
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A Renormalized excitonic method for an 
(A - B) n chain 

One has two types of blocks A and B. The ground state 
and lowest excited eigenfunctions for each block are given 
by 



Ha\^a) 
HbWb) 



= ^bWb), 
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The ground state of the chain is represented by 

*° = IWII*V 

» 3 

Its energy will be 

(*°\H e "\*°) = n{E A + E% ) + n(« AB + 

The interaction energies between two adjacent blocks is 
given by the knowledge of the exact energies of the AB 
and BA dimers 

H AB \* AB ) = E° AB \* AB ), E% = E A + E° B + v AB . 

For the description of excited states one needs to estimate 
the effective interaction between a local excited state and 
the neighbor ground states and the integral responsible for 
the transfer of excitation. These informations are obtained 
from the excited solutions of the dimers. One should es- 
pecially consider the two eigenstates 

Hab^Xb) = Eab\*ab), 

h ab \^; b ) = e*; b \^; b ), 

having the largest projections I&Xb) an d I^ab); on the 
model space spanned by ip A ip B and ip A ip B , which can be 
written after orthogonalization as 

\n B ) = m^D + M°Ar B ), 

\Hb) = -/#aV&> + A|^>- 
It results that 

(r A K\H eff \r A ^B) = \ 2 e ab +v 2 e*; b 

= E* A + E° B +v {A * )B , 

(r A r B \H eff \r A r B ) = v 2 e* ab + \ 2 e*; b 

= E A + E* B + u J 4(s*) ; 

(rA^ B \H eff \i> A r B ) = (e* ab - e*; b )\» = h AB . 

For the periodic system the delocalized excited states will 
be respresented as linear combinations of locally excited 
states on either A or B blocks 

jjtn i 

The energies of and W Bn are given by 

(& A jH eff \&XJ - (&°\H eff \&°) = 
E * A - E° A +v {A , )B - v AB , 
(& Bn \H eff \& Bn ) - (&°\H eff \&°) = 
E* B - E% + - v BA , 

This locally excited state are coupled with the states lo- 
cally excited on the adjacent B blocks 

(njH ef, \nj = hAB. 



For the lowest state of the lattice, corresponding to k =0, 
the delocalized excited states can be written as a linear 
combination of 



and 



solution of a 2 x 2 matrix whose elements are 



(Kh =0 # e// (KhJ 



H eff 



E* A — E A + 2^^.)^ - V AB ^, 



H 



eff 



H 



eff 



(*f)t=0y 

E* B - E B + 2{Vi yB ,- )A - V BA j , 



k =0 



H 



eff 
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